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O ; We reinvestigate a simple relation between the semileptonic CP asymmetry a s sl , the decay 
rate difference Ar s , the mass difference AM S and extracted from the angular analysis 
^ ■ of the decay B s — > ip<p, which is regularly used in the literature. We find that this relation 
! is not suited to eliminate the theory prediction for Ti 2 , it can, however be used to determine 
the size of the penguin contributions to the decay B s — > ipcf). Moreover we comment on the 
current precision of the theory prediction for r 12 . 
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1 Introduction 



Currently we have some hints for deviations of experiments from standard model predictions 
at the three sigma level both in the B&- and the S s -mixing system [U EJ EJ HI EJ Ej . In the 
-Bs-system some of the central values differ largely from the corresponding standard model 
values - here LHCb will clearly tell us till the end of 2011 whether these large deviations are 
realized in nature. In view of the expected precision of the coming data, it is mandatory to 
try to achieve the same precision in the theory predictions. In this paper we revisit a simple 
relation between four mixing observables and show that this relation can be badly violated, 
due to neglecting some small quantities. 

After reviewing the mixing formalism in Section 2.1, we show how new physics affects different 
phases arising in the mixing of neutral B-mesons. In Section 2.3 we derive the relation and we 
test it in the following section. In Section 3 we discuss the accuracy of the theory prediction 
for Ti2. Finally we conclude in Section 4. 



2 The Relation 

2.1 Mixing formalism 

Mixing of neutral B mesons is described by the off-diagonal elements Mf 2 and r^ 2 of the mass 
and decay rate matrix. The following mixing observables are determined in experiments: the 
mass difference 

AM, := M%-Ml = 2\M\ 2 \ (l + sin 2 <J> q + ...), (2.1) 



the decay width difference 



: r| - r% = 2|r? 2 | cos0, ( i - \^- 2 sin 2 <j> q + ... ; 



8|M« 2| 



and the flavour specific asymmetry 



if l / r q 

sin (6.) +0 — 
\M 12 \ lVV V^12 



' sin +° ( TFT ) » (2-4) 



(2.5) 



where <p q = arg(— Mf 2 /rf 2 ). A typical example of a flavor-specific decay is a semileptonic 
decay, therefore this asymmetry is also called semileptonic CP-asymmetry, a q sl = a q j s . In the 
standard model rf 2 /Mf 2 is of the order of 5 • 10~ 3 , so one can safely neglect {T\ 2 / Ml 2 ) 2 ■ The 
standard model values of these quantities were recently updated in [7] using the NLO-QCD 
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calculations from 019]. 

AMf = (17.3 ± 2.6) ps~\ AM| M = (0.54 ± 0.09) ps" 1 , (2.6) 

— = 0.137 ±0.027, — = (4.2 ± 0.8) • 10" 3 , (2.7) 

a ^ SM = (1.9 ± 0.3) • 10~ 5 , a d f f M = - (4.1 ± 0.6) • 10" 4 , (2.8) 
0SM = 0.22° ±0.06°, = -4.3° ±1.4°, (2.9) 

A b f M = 0.506aj SM + 0.494a*f M = (-2.0 ± 0.3) • 10" 4 , (2.10) 

<f M ~ 4 SM = ( 4 -3 ± 0.7) ■ 10" 4 . (2.11) 
All these values were obtained with the input parameters taken from [3]. 

2.2 New Physics contributions to mixing 

In the general case of new physics being present in S-mixing we can write model-indepently 

M\ 2 = Mlf M ■ A q = Mlf M ■ \A q \e 1 ^ , (2.12) 

r? 2 = r^ SM -A g = r^ SM -|A g | e -^ , a,»i. (2.13) 

Then the mixing phase can be decomposed as 

0, = 0f + 0? + 0? • (2-14) 

The standard model part 0^ is tiny in the case of B s mesons (Eq. (12. 9p ). the new physics 
contribution to Mf 2 is denoted by 0^ and the hypothetical new physics contribution to 
r^ 2 , which is strongly constrained by different well-measured observables is denoted by 4>q- 
Because of these contraints we will neglect new physics contributions to Y\ 2 in the following^] 
A related quantity arises in the angular analysis of the decay B s — > ifxj), which is sometimes 
confused with cj) s [10]. 2(3 S := — aig[(V t bV t *) 2 / (V c bV* s ) 2 } is the phase which appears in b — > ccs 
decays of neutral B-mesons taking possible mixing into account, so e.g. in the case B s — > 
iy t bV* s ) 2 comes from the mixing (due to M i2 ) and (V c bV* s ) 2 comes from the ratio of b — > ccs 
and b — > ccs amplitudes. Sometimes (3 S is approximated (using the PDG convention for the 
CKM elements!) as 2(3 S pa — arg[(y t bV t * s ) 2 ] ~ — arg[(V t *) 2 ] - the error due to this approximation 
is on the per mille level. The standard model value for this angle reads [3] 

2(3 S = (2.1 ±0.1)° . (2.15) 

As mentioned above <fi s : = arg[— M 12 /r 12 ] is the phase that appears e.g. in a s j- s . In M 12 
we have again the CKM elements {V t bV t * s ) 2 , while we have a linear combination of {V c bV* s ) 2 , 
VcbV* s VubVus an d (VubVus) 2 i n ^12- Neglecting the latter two contributions - which is not 



1 Taking a non-negligible value of 0^ into account, will only strengthen our arguments. 



2 



justified - would yield the phase 2/3 s . Numerically 0f is one order of magnitude smaller 
than 2(3 S . New physics contributions to mixing, i.e. to M{ 2 alters (3 S as 

- 2(3 S -)■ (f>f - 2(3 S . (2.16) 

Taking into account possible penguin contributions both in the standard model and beyond 
one gets 

- 2/3, + §r ng ' SM -> <f>f - 2(3 S + 6* ens ' SM + 5f° ng ' NP . (2.17) 

The penguin contributions are typically expected to be small [11]. An interesting exception 
of the decomposition in Eq.f l2.17]) is given by the standard model with four generations. In 
this particular model a sizeable deviation of — 2f3 s from its standard model value is possible, 
see e.g. [T2l [13] for some bounds on the CKM-elements in this model. 
Sometimes in the literature (e.g. [HI [151 E]) the following quantity is used 

= sin (-<f>? + 2ft - 5f ng ' SM - 5f ng > NP ) ■ (2. 18) 

A model independent fit [3] for new physics in 5-mixing, gives the following currently allowed 
range [7] 

SW = q - 78 -oaI > (2-19) 
which has to be compared with the SM value 

Sim = 0.036 ± 0.002 . (2.20) 

Penguin contributions have been neglected in Eq.f l2.19l) and Eq. (12.20]) . In the literature it is 
sometimes argued, that if the new physics contribution is sizeable, then we can approximate 

4> S s M + 4>? + 4>?~4>?, (2-21) 

0f - 2(3 S + 5^ cns ' SM + 5 p s cne ' NP « </>? , (2.22) 

since the standard model phases and the possible penguin contributions are very small. 



2.3 Deriving the relation 



This approximation (Eq.f l2.21]) and Eq. fl2.22p ) was used e.g. in [TU [151 ES1 EE] to derive a 
simple model independent relation between observables in the mixing system. 



'si 



1 12 



sin 



(0s SM + <f>t) 



M{ 2 

2|ry CO s {<p s s M + <p?) 



2 Mi 



tan 



121 
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with the correction factor 5 

5 (/3 S SM , 0f M , 5f ng ' SM + 6? ng ' NP , 0f ) 



tan 



lSM 



(2.23) 



tan (-2/3fM + 0A + ^.sm + ^.np 



.(2.24) 



Since recent fits [3j 0] of the new physics contribution to the B s mixing phase give relatively 
large values 



0? = (-52±!) 



)° at 95 % C.L. [3], 



(2.25) 



(-40 ± 16)° at 68 % C.L. @], 



(2.26) 



,s 



it seems to be justified to neglect all other contributions -2/3f M and < jpeng,SM + ( jpeng,NP) 

in Eq. fl2.24p and to use 5=1. This was done often in the literature and one gets 



This relation corresponds to Eq.(23) in [13], to Eq.(A.l) in [15], to Eq.(55) in [16] and to 
Eq.(3) in [T7]. Several years earlier similar relations to Eq. fl2.27p were derived e.g. in [TB~| 
[1911201 [21]. In Eq.(7.10) of [H] and Eq.(10) of [19] also the same approximations (Eq.fl22D 
and Eq.( l2.22p ) were used, but in the final formula there was still theory input on r 12 /M 12 . 
Therefore one expects also a correction factor 5, which can be obtained from 5 by replacing 
the tangent with the sine. Therefore 5 deviates less from one than 5. For Eq.(26) of [20J only 
the approximation in Eq. (12.21 p is used, which is expected to work well. For Eq.(3.53) of [21] 
none of the above approximations was used. 

Similar relations without the approximations in Eq. (l2.2ip and Eq. (12.22]) were presented e.g. 
in Eq.(10) of [H], Eq.(3.5) in [22] and in [TB]. Kagan and Sokoloff [IB] also discuss deviations 
from the relation Eq. (12.27]) by expanding to first order in the penguin contributions and the 
small phases, see Eq. (110-113) of (16"] . 

2.4 Testing the Relation 

In order to test the approximations made for deriving Eq. (l2.27p we plot the correction 
factor <5(/3 s SM ,0^ M ,5P cn g' SM + (5P eng ' NP ,0f) as a function of <pf for four different values of 
the unknown penguin contributions #P en g' SM + ^peng.NP _ qo^ _2° j —5° and —10°. Compar- 
ing with [11] the second and the third values seem to be conservative and realizeable in 
nature^). <5P en S' SM + ,5P cn g: NP — 1 o° corresponds to large penguin contributions (see how- 
ever [23], where such a possibility is even not excluded within the standard model), while 
£peng,SM _|_ ^peng.NP _ go correS p 0n( j s ^ Q the hypothetical case that standard model penguins 
and new physics penguins cancel exactly. /3f M and 0f M are fixed to their standard model 
values given Eq. (l2.15p and Eq.( )2.9p . 

2 More precise data might shrink further the allowed range for penguin contributions. 




(2.27) 
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As expected the relation in Eg. (12.271) does not work at all for small new physics contributions 
(p^, i.e. 0^ ~ or (fif « it, but one can see from the plot, that the relation in Eg. (12.271) does 
also not work for large new physics contributions. 

To investigate the accuracy of Eg. (l2.27p further let us first discuss the case of a very large 
value of the new physics mixing phase (pf, where the approximations from Eg. (12. 21 p and 
Eg. (l2.22p are expected to work best. Taking the large value from Eg. (12 .25]) one gets the 
following values for 5: 





£pcng,SM+iNP _ qo 


^pcng,SM+JNF _ 


£pcng,SM+NP _ _go 


£pcng,SM+NP _ _1Q° 


5 


1.08638 


1.16571 


1.29463 


1.54297 



If there would be no penguin pollution at all, the relation Eg. (l2.27p would have a maxi- 
mal accuracy of about 10%, even for moderate penguin pollutions Eg. (l2.27p is violated by 
20% — 30% (for smaller values of (pf, the situation is much worse), while for large penguin 
contributions we have a violation of Eg. (l2.27p of more than 50%. So even in the ideal case 
of a very large value of the new physics mixing phase 4>f a violation of up to (9(50%) can 
easily occur. If 0^ does not have the value from Eg. (l2.25p . then the relation in Eq. fl2.27p can 
of course be violated much stronger than (9(50%). 

Next we show a numerical example, to illustrate that a measurement of a deviation of 5 from 
one has to combined with theoretical input in order to obtain information about the phases. 
Let us assume that all the quantities in Eg. (12.271) are measured and e.g. 5 = 1.6 will be found 
experimentally. Such a value can correspond to very different values of new physics phases 





£peng,SM+NP _ qo 


£peng,SM+NP _ 


£pcng,SM+NP _ _go 


£peng,SM+NP _ _^Qo 


<t>f 


0.10457 


0.200114 


0.352488 


0.69367 



Without including further input (like the theoretical determination of T 12 ) the new physics 
phase can lie somewhere between 0.1 und 0.7 and we cannot distinguish new physics effects 
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in mixing from QCD effects like penguin contributions. 

T or 0f 



Including however theoretical input like Tf^ or <fr^ M we can determine all phases, e.g. 



• a s sl from experiment and Tflf and from theory allows us to extract (pf. This can 
then be combined with the experimental value of or 5 to obtain — 2f3 s + $P en g> SM + 
£Peng,NP_ y[ os i np models (except e.g. the SM4) do not change the value of — 2j3 s and 
therefore we can use the SM prediction for — 2f3 s to get a precise determination of the 
penguin contributions 

• The experimental values of and 5 can be used to extract 0™ + (fif . With the 
theoretical value of 0f M the new physics phase (f>f can be extracted. As above one can 
now extract the penguin contributions. 

To summarize the results of this section: The relation Eq. ( 12.271) can be strongly violated 
due to the small difference between and — 2f3 s and non-vanishing penguin contributions. 
Therefore a violation of Eq. fl2.27p . i.e. 8^1, does not provide a test of our theoretical 
framework for deriving T 12 and M 12 . Eq. (12.271) can not be used to eliminate the theory 
prediction for T 12 . The exact ralation Eq. (l2.23p (with the correction factor S) can however 
be used to extract the desired information: Using the theory prediction for Ti2 one can extract 
cj)f and the size of the penguin contributions ^P eng ' SM + 5P eng,NP , which is an important result. 

3 Comment on the theoretical accuracy of 

The DO collaboration measured [24"l [25] a very large value for the Dimuonasymmetry 

A b,SM = _o.o0957 ± 0.00251 (stat) ± 0.00146(syst) , (3.28) 

which differs 3.2 a from the standard model prediction [3, [9], or Eq.( 12.10l) . This result 
triggered a lot of theoretical interest, see e.g. [12], [26] (due to a lack of space we quoted 
only papers from the first two months after the DO result appeared). Allowing only for 
new physics in Mf 2 (Eq. fl2.12l) ) one gets the following relation for flavor-specific/semileptonic 
CP-asymmetries 

|r? 9 1 sin(0;>M + 0A) 

a q fs = 1 % V , 9 A , q> . (3.29) 

fs \Mfi M \ \A q \ 1 ; 

From this relation one can derive a bound on the maximal value of the dimuon-asymmetry 

^ MAX ~(-5±1)-10~ 3 , (3.30) 

which is about 1.5 a below the experimental value in Eq. fl3.28p . 

Due to this discrepancy (although the statistical significance is only 1.5 a), it was suggested 
[27] that new physics might also act in T\ 2 , c.f. Eg. (12.131) or that the theory prediction for 
r\ 2 might be affected by non-perturbative effects. One possibility to circumvent hypthetical 
problems with T\ 2 would be the elimination of the corresponding theory prediction with the 
help of Eq.f l2.27p . as suggested e.g. in [17]. But as explained above Eq. fl2.27p can not be used 
without theory information on T 9 l2 . 

In order to shed some light on the necessity of the elimination of the theory prediction for 
Ti2 we review here its theory status. T\ 2 has three contributions 

^12 = ^c^cc + 2A c A u r ttc + X U T UU . (3.31) 
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T xy corresponds to a box diagram with internal x- and y— quarks and Aj = V* s Vib. We 
investigate now the expected expansion parameter in the Heavy Quark Expansion for the 
individual Y xy . 

T cc : This contribution dominates by far |r 12 | and -Re(T 12 ) and therefore describes Ar s . 
Since we have now two charm quarks in the intermediate state, the expansion parameter 
of the Heavy Quark Expansion is not the inverse bottom mass but a reduced mass: 

A A A 1 , 
> , = ; • 3.32 

Using pole masses for the quarks one gets an expansion parameter of about 1.3A/m&. 
It is however well-known that the use of the pole mass suffers from considerable uncer- 
tainties related to renormalons. Using instead the method and parameters (MS-values 
at the same scale for the quark masses in z, which corresponds to summing up loga- 
rithms of the form z In z to all orders) which were used in [7] we get as an expansion 
parameter of the HQE 

A 1 A , 

^1.11 ±0.01) — , (3.33) 



m b a/1 - Az m b 

that is almost identical to A/m&. So this simple dimensional estimate does not indicate 
any problems concerning the convergence of the HQE. 

Moreover the validity of the HQE for T cc can be tested directly by comparing theory 
and experiment for Ar s and indirectly by the lifetime ratio T^ s /r^ d , because a very 
similar contribution arises in the theoretical determination of the lifetime of the B s 
meson, see e.g. [28]. Currently no deviation from the standard model predictions are 
seen [7j, but more precise experimental numbers for Ar s and tb s are very desireable. 

F U c,uu '■ These two contributions give the dominant contribution to Jm(r 12 ) and are therefore 
important for a s i. Since we now have at most one charm quark and else only light 
up-quarks as internal quarks, naive power counting shows that the HQE is given as an 
expansion in the inverse heavy b-quark mass, which is expected to be well-behaved. 
A very similar contribution arises in the theoretical determination of the lifetime of 
the Bd and B + mesons. Theory and experiment agree well for the ratio r B +/rB d [7J, 
although the theoretical precision is strongly limited by a lack of knowledge of the 
arising non-perturbative parameters. 

As we have shown, dimensional estimates do not indicate a breakdown of the convergence 
of the HQE, but instead of dimensional estimates it is much more instructive to determine 
explicitly the size of all the corrections to T 12 . We can write 

Ar s = Ar° (1 + 5 Latticc + 5 QCD + 5 HQE ) . (3.34) 

Ar° is the theory prediction in LO-QCD, LO-HQE (i.e. only contributions of dimension 6) 
and with all bag parameters set to one. S Lattlce corresponds to the deviation of the lattice 
results for the bag parameters from one, <5^ CD corresponds to the NLO-QCD corrections and 
^hqe ^ Q highgr orders in the HQE. With the numerical values used in [7] we get the 
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following results El 



0.142ps 



-i 



(3.35) 
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Lattice 



0.14 



(3.36) 



§QCD 



0.06 , 



(3.37) 



£HQE 



0.19 . 



(3.38) 



All corrections are negative and smaller than 20% = 1/5. So the direct calculation of the 
first order corrections suggests that the HQE is well-behaved. This can be compared with 
the status of 2004 [29] , where considerable larger uncertainties were still present in the theory 
prediction for Ti 2 . 

4 Conclusion 

We have investigated the relation 



and shown that it receives verly large corrections in dependence of the value of the new 
physics phase in mixing. Even in the ideal case of very large values of the relation 
can be violated up to 50 %. Therefore it can not be used to eliminate the theory prediction 
for r^ 2 , but including the correction factor 5 from Eq. (12.24] ) it can be used to determine the 
size of the penguin contribution to the decay B s — > ip<p, which is a very important task. A 
sizeable penguin contribution can also be a signal for new physics. 

We also have reinvestigated the accuracy of the theory determination of Ti2 and found no sign 
for unexpectedly large corrections within the framework of the HQE. Comparision between 
experiment and theory predictions within the framework of the HQE shows a good agreement, 
with one exception: the central value of the dimuon-asymmetry measured by DO is 1.5 a above 
the theory bound. Although this discrepancy is statistically not significant, more precise data 
for the dimuon asymmetry from TeVatron would be very helpful. Moreover with the expected 
new data on AT S and Tb s - in particular from LHCb - soon much more profound conclusions 
about the value of Ti2 can be drawn. 
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